First Principles Calculation of Anomalous Hall Conductivity in Ferromagnetic bcc Fe 
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We perform a first principles calculation of the anomalous Hall effect in ferromagnetic bcc Fe. 
Our theory identifies an intrinsic contribution to the anomalous Hall conductivity and relates it to 
the k-space Berry phase of occupied Bloch states. The theory is able to account for both dc and 
magneto-optical Hall conductivities with no adjustable parameters. 
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In many ferromagnets the Hall resistivity, pu, exhibits 
an anomalous contribution proportional to the magne- 
tization of the material, in addition to the ordinary 
contribution proportional to the applied magnetic field, 
p„ = RoB + RAttM [11113. The anomalous Hall effect 
(AHE) has played an important role in the investigation 
and characterization of itinerant electron ferromagnets 
because Rs is usually at least one order of magnitude 
larger than the ordinary constant Rq. Although the effect 
has been recognized for more than a century |3] it is still 
somewhat poorly understood, a circumstance reflected 
by the controversial and sometimes confusing literature 
on the subject. Previous theoretical work has failed to 
explain the magnitude of the observed effect even in well 
understood materials like Fe Ql ■ 

Karplus and Luttinger [sl pioneered the theoretical in- 
vestigation of this effect, by showing how spin-orbit cou- 
pling in Bloch bands can give rise to an anomalous Hall 
conductivity (AHC) in ferromagnetic crystals. Their con- 
clusion was questioned by Smit 0], who argued that Rg 
must vanish in a periodic lattice. Smit proposed an al- 
ternative mechanism, skew scattering, in which spin-orbit 
coupling causes spin polarized electrons to be scattered 
preferentially to one side by impurities. The skew scat- 
tering mechanism predicts an anomalous Hall resistivity 
linearly proportional to the longitudinal resistivity; this 
is in accord with experiment in some cases, but an ap- 
proximately quadratic proportionality is more common. 
Later, Berger proposed yet another mechanism, the 
side jump, in which the trajectories of scattered elec- 
trons shift to one side at impurity sites because of spin- 
orbit coupling. The side jump mechanism does predict 
a quadratic dependence of the AHC on the longitudi- 
nal resistivity. However, because of inevitable difficulties 
in modeling impurity scattering in real materials, it has 
not been possible to compare quantitatively with experi- 
ment. It appears to us that the AHE has generally been 
regarded as an extrinsic effect due solely to impurity scat- 
tering, even though this notion has never been critically 



tested, and that the intrinsic contribution initially pro- 
posed by Karplus and Luttinger has been discounted. 

Several years ago, the scattering free contribution of 
Karplus and Luttinger was rederived in a semiclassical 
framework of wavepacket motion in Bloch bands by tak- 
ing into account Berry phase effects 8, 9] . According to 
this work, the AHC in the scattering free limit is a sum of 
Berry curvatures (see eqs.© and © below) of the occu- 
pied Bloch states lOj. Recently, Jungwirth et al. [llLll^ 
applied this picture of the AHE to (III,Mn)V ferromag- 
netic semiconductors and found very good agreement 
with experiment. (III,Mn)V ferromagnets are unusual, 
however, because they are strongly disordered and have 
extremely strong spin-orbit interactions. In this Letter, 
we report on an evaluation of the intrinsic AHC in a 
classic transition metal ferromagnet, bcc Fe. Our calcu- 
lation is based on spin-density functional theory and the 
LAPW method. The close agreement between theory 
and experiment that we find leads us to conclude that 
the AHC in transition metal ferromagnets is intrinsic in 
origin, except possibly at low-temperature in highly con- 
ductive samples. 

We begin our discussion by briefly reviewing the semi- 
classical transport theory. By including the Berry phase 
correction to the group velocity 0, 0| , one can derive the 
following equations of motion: 
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is the Berry curvature of the Bloch state de- 
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with Unk being the periodic part of the Bloch wave in the 
?ith band. We will be interested in the case of B = 0, 
for which £„(k) is just the band energy. The distribu- 
tion function satisfies the Boltzmann equation with the 
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usual drift and scattering terms, and can be written as 
/„(k) + (5/„(k), where /„ is the equiUbrium Fermi-Dirac 
distribution function and (5/„ is a shift proportional to 
the electric field and relaxation time. The electric cur- 
rent is given by the average of the velocity over the dis- 
tribution function, yielding to first order in the electric 
field [13 

n n 

(3) 

The same expressions can be derived from the Kubo 
linear-response-theory formula for the conductivity ma- 
trix. The first term is the anomalous Hall current orig- 
inally derived by Karplus and Luttinger, but never pre- 
viously evaluated. In the second term, apart from the 
longitudinal current, there can also be a Hall current in 
the presence of skew scattering because the distribution 
function can acquire a shift in the transverse direction. 
This skew scattering contribution to the Hall conductiv- 
ity should be much smaller than, but crudely propor- 
tional to, the longitudinal conductivity and can be iden- 
tified, when it is dominant, by the traditional test, i.e. 
Pxy oc pxx- It will have a larger relative importance when 
Uxx is large, i.e. in pure-crystals at low temperatures. 

We now discuss our scheme for calculating the Berry 
curvature and the AHC. For a cubic material with mag- 
netization aligned along [001], only the ^-component 
i7^(k) ^ 0. In our calculation, we find it convenient 
to use a different but equivalent expression for the Berry 
curvature that arises naturally from the Kubo formula 
derivation 141, 
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where En = fi^m and w's are velocity operators. In 
the relativistic formulation, the ^/'nk are four-component 

Bloch wave functions, and w — ^ ), with a being 

the Pauli matrix and c the speed of light. It will also 
be instructive to introduce the sum (for each k) of Berry 
curvatures over the occupied bands: 



0^(k)=^/„f]J^(k). 
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Then the intrinsic AHC is an integration over the Bril- 
louin zone (BZ): 
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The recent development of highly accurate ab initio 
electronic structure calculation methods enables us to 
complete the work of Karplus and Luttinger by eval- 
uating their intrinsic Hall conductivity and comparing 



it with experiment. We employ the full-potential lin- 
earized augmented plane- wave method [l5| with the gen- 
eralized gradient approximation (GGA) for the exchange- 
correlation potential Fully relativistic band cal- 
culations were performed using the program package 
WIEN2K A converged ground state with magneti- 
zation in the [001] direction was obtained using 20,000 
k points in the first Brillouin zone and K^^xRut = 
10, where Rmt represents the muffin-tin radius and 
^max the maximum size of the reciprocal-lattice vectors. 
In this calculation, wavefunctions and potentials inside 
the atomic sphere are expanded in spherical harmonics 
Yim{0, 4>) up to ^ = 10 and 4, respectively, and 3s and 3p 
semi-core local orbitals are included in the basis set. The 
calculations were performed using the experimental lat- 
tice constant of 2.87 A. The spin magneton number was 
found to be 2.226, compared to the experimental value 
of 2.12 as deduced from measurements of the magnetiza- 
tion Jl^ and of the g (= 2.09) factors. The calculated 
energy bands are shown in Fig.l, and are very similar to 
those obtained in Ref. |0| . If the spin-orbit interaction 
is parameterized as C 1 ' its strength ^ is found to be 
approximately 5.1 mRy from the band splitting near the 
H point and the Fermi energy. 
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FIG. 1: Band structure of iron along high symmetry lines 
in the Brillouin zone. The magnetization direction is along 
[001]. 

After obtaining the self-consistent potential with 
20,000 k points, we calculated the Berry curvature with 
several larger sets of k-points in order to achieve the con- 
vergence for Gxy shown in Fig. 2. The Monkhorst-Park 
special-point method |20j was used for the integration in 
Eq.®. To go beyond 2,000,000 points, we adopted a 
method of adaptive mesh refinement, i.e., when r2^(k) is 
large at a certain k point, we construct a finer mesh by 
adding 26 additional points around it. This procedure 
yields a converged value of Gxy = 751 (r2cm)~^ at zero 
temperature (using a step function for the Fermi-Dirac 
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FIG. 2: Anomalous Hall effect vs. 5 with different numbers 
of k points in full Brillouin zone. Here S is introduced by 
adding 5^ to the denominator in Eq.Q- The solid lines were 
obtained by an adaptive mesh refinement method. 

distribution) and a slightly smaller value of a^y = 734 
(f2cm)~^ at room temperature (300 K). Our result is in 
fair agreement with the value 1032 (i}cm)~^ extracted 
from Dheer's data on iron whiskers |2l| at room temper- 
ature. 

The slow convergence is caused by the appearance of 
large contributions to of opposite sign which occur 
in very small regions of k-space. Spin-orbit effects are 
small except when they mix states that would otherwise 
be degenerate or nearly degenerate, and even then, those 
mixed states will contribute nearly canceling contribu- 
tions to . Only when the Fermi surface lies in a spin- 
orbit induced gap is there a large contribution. This can 
be seen in Fig. 3 where the Berry curvature along lines in 
k-space is compared with energy bands near Ep and in 
Fig. 4 where it is compared with the intersection of the 
Fermi surface with the central (010) plane in the Brillouin 
zone. 

In order to further understand the role of spin-orbit 
coupling in the AHE, we artificially varied the speed of 
light, thereby changing the spin-orbit coupling strength 
^ oc c~^. As shown in Fig. 5, axy is linear in ^ for small 
coupling, but not for large coupling. For iron, nonlin- 
earities become significant for ^/£^q > 1/2, which means 
that the spin-orbit interaction in iron cannot be treated 
perturbatively. 

So far we have discussed only the dc-AHE. It is 
straightforward to extend our calculation to the ac Hall 
case by using the Kubo formula approach: 

<j{uj)^y = ^ / (-^"^k - /n',k) (7) 



FIG. 3: Band structure near Fermi surface (upper figure) and 
Berry curvature Q,^ (k) (lower figure) along symmetry lines. 
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FIG. 4: (010) plane Fermi-surface (solid lines) and Berry 
curvature — r2^(k) (color map). —Q.z is in atomic units. 



where e is a positive infinitesimal. In the upper panel in 
Fig. 6, we show results for the imaginary part of oJOxy as 
a function of frequency that are in agreement with earlier 
calculations [2^. Experimental results are in excel- 
lent agreement below 1.7 eV but become smaller at higher 
energies. In the lower panel of the figure, the real part of 
the Hall conductivity, obtained from the imaginary part 
by a Kramers-Kronig relation, is shown as a function 
of frequency. The dc limit result, (j{u} = Q)xy = 750.8 
(ficm)"^, is essentially identical to that obtained from 
Eq.(6). Despite the small discrepanc y w ith theory in the 
dc limit, the experimental point • |21| seems to agree 
rather well with the overall trend of the frequency de- 
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FIG. 5: Calculated anomalous Hall conductivity (open cir- 
cles) vs. spin-orbit coupling 5, is spin-orbit coupling 
strength of iron. The line is a guide to the eye. 



pendence of the calculated AHC. 
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FIG. 6: Frequency dependence of the Hall conductivity at 
zero temperature. In the upper panel, the calculated imagi- 
nary part of ujaxv (solid curve) is compared with experimental 
results o Ref. (|2^). In the lower panel, the real part of Oxy 
is shown together with the dc experiment value • extracted 
from Ref ('2TI). 

In conclusion we have shown that the AHC of bcc 
Fe, and presumably all other transition metal ferromag- 
nets, is primarily intrinsic. (The only previous evalua- 
tion of the AHC of which we are aware found that 
axy — 20.9(r2cm)^^) for Fe. The remaining discrepancy 
between theory and experiment is likely due to short- 
comings of the GGA, neglect of scattering effects, and 
experimental uncertainties. 
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